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ABSTRACT 

The performance of three recently proposed second-order closure models is tested in 
benchmark turbulent shear flows. Both homogeneous shear flow and the log-layer of an 
equilibrium turbulent boundary layer are considered for this purpose. An objective analysis 
of the results leads to an assessment of these models that stands in contrast to that recently 
published by other authors. A variety of pitfalls in the formulation and testing of second- 
order closure models are uncovered by this analysis. 
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INTRODUCTION 


The need for advanced turbulence models to reliably compute the complex aerodynamic 
flows of technological interest has led to a resurgence of interest in second-order closure mod- 
els. Consequently, the recent papers by Shih, Chen and Lumley 1 and Shih and Lumley , 
which reported tests of second-order closure models in turbulent shear flows, attracted our 
attention. In these papers, results were presented which appear to indicate that the Shih- 
Lumley model 3 performs better than other recently proposed second-order closures in homo- 
geneous shear flow as well as in other boundary-free turbulent shear flows. However, our own 
comparative studies of second-order closure models in benchmark turbulent shear flows have 
yielded a different picture. The purpose of the current paper is to present these alternative 
results for comparison. 

The predictions of three second-order closure models recently proposed by Shih and 
Lumley 2 ’ 3 , Fu, Launder and Tselepidakis 4 and Speziale, Sarkar and Gatski 5 will be com- 
pared in two benchmark turbulent flows: homogeneous shear flow and the log-layer of an 
equilibrium turbulent boundary layer. These flows are selected since the former constitutes 
a basic building-block free turbulent shear flow whereas the latter serves as a cornerstone for 
the calculation of practical wall-bounded turbulent flows of engineering interest. Particular 
attention will be paid to evaluating the ability of each model to accurately predict the equi- 
librium values for the Reynolds stress anisotropies. However, for the case of homogeneous 
shear flow, model predictions for time evolving fields will also be compared. Objective means 
for evaluating the performance of the models will be provided and pitfalls in the formulation 
and evaluation of models are uncovered that have led to previously published assessments 
that are misleading. 

THE TURBULENT SHEAR FLOWS TO BE CONSIDERED 

We will consider incompressible turbulent shear flows with the mean velocity gradient 
tensor 

P- = SSnSfi ( 1 ) 

OXj 

where Sij is the Kronecker delta and S is shear rate. In homogeneous shear flow, the shear rate 
S is constant and is applied in an unbounded flow domain yielding spatially homogeneous 
turbulence statistics. For this, as well as any homogeneous turbulence, the Reynolds stress 
tensor r,j is a solution of the transport equation 6 



where T,y = u,Uj and 


( dui duj\ 

' j ~ p {d^ + d^)' 


duj duj 
E ' 3 V dx k dx k 
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are, respectively, the pressure-strain correlation and dissipation rate tensor (here, p is the 
fluctuating pressure, u t is the fluctuating velocity, v is the kinematic viscosity, and an overbar 
represents an ensemble mean). Since 


£,j — n ij 


( 4 ) 


where II, y = $,y — ij£,j given that p £ ij Is the deviatoric part of the dissipation rate tensor 
(e = closure is achieved once models for II, y and e are provided. In most existing 

second-order closure models, £>c,y is neglected while II, y and e are modeled in the general 
form: 

dv k 
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dv, 
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are the turbulent kinetic energy and Heynolds stress anisotropy tensor, respectively. Here, 
Aij and Mijki are dimensionless tensor functions of and possibly the turbulence Reynolds 
number R t = K 2 jue\ C c \ and C C 2 are either constants or functions of the second and third 
invariants (II , III) of as well as R t . The full form of the three models to be considered 
- the Shih-Lumley (SL) model, the Fu, Launder and Tselepidakis (FLT) model and the 
Speziale, Sarkar and Gatski (SSG) model - are provided in the Appendix. 

For homogeneous shear flow, each of the models - consistent with physical and numerical 
experiments - predict that the anisotropy tensor and shear parameter SKje achieve equi- 
librium values that are independent of the initial conditions (see Tavoularis and Corrsin 7 , 
Rogers, Moin and Reynolds 8 and Tavoularis and Karnik 9 ). This equilibrium state is associ- 
ated with solutions where = 0 or equivalently, 


rij = {V- ef-g 


( 8 ) 


where V = —Tijdvi/dxj is the turbulence production. The substitution of (1) and (8) into 
(2) yields the system of algebraic equations (see Abid and Speziale 10 ) 


min 
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where 11*^ = Tl,j/ SK is, for the turbulent shear flows to be considered, a function of TV IK 
and V je whose specific form depends on the pressure-strain model chosen. In deriving (9), 
the identity 


( 10 ) 


V _ _n 2 / SK\ 
e ~ K V e ) 

has been used. Once the equilibrium value of V je is specified (and Eq. (7) is utilized), it 
is straightforward to obtain the equilibrium values of from a numerical solution of (9). 
These equilibrium values are determined exclusively by the pressure-strain model. 

In the logarithmic region of an equilibrium turbulent boundary layer, the mean velocity 
gradient tensor is of the general form (1) and there is a production-equals-dissipation equi- 
librium where the turbulent diffusion terms vanish in the Reynolds stress transport equation. 
Consequently, (9) yields the equilibrium Reynolds stress anisotropies for the log-layer when 
V je — 1 (see Abid and Speziale 10 ). In contrast to this, homogeneous shear flow achieves 
an equilibrium state where V/e > l (physical and numerical experiments 7 ” 9 have indicated 
that 1.4 < V/e < 1.8). 


DISCUSSION OF RESULTS 

A comparison of the results predicted by the pressure-strain models of Shih and Lumley 2,3 , 
Fu, Launder and Tselepidakis 4 and Speziale, Sarkar and Gatski 5 will now be made. In 
Figures 1(a) and 1(b), the model predictions for the norm of the slow and rapid parts of the 
pressure-strain correlation are compared with the DNS results of Rogers, Moin and Reynolds 8 
for homogeneous shear flow (Hy = II + 11,-j^ where the superscripts (5) and ( R ) denote, 
respectively, the slow and rapid parts; (n^II.j) 1 / 2 denotes the L 2 norm of Ily)). These 
results are similar to those used by Shih and Lumley 2 to conclude that the Shih-Lumley 
model performs the best of these three models and that the SSG model performs poorly. 
Such a conclusion is highly misleading. The SSG model is a model for the total pressure- 
strain correlation and not for its separate slow and rapid parts. The standard hierarchy of 
pressure-strain models (5) are only theoretically justified for homogeneous turbulent flows 
that are near equilibrium (as shown by Speziale, Gatski and Sarkar 11 , both the slow and rapid 
parts of the pressure-strain correlation depend nonlinearly on the mean velocity gradients at 
retarded times in non-equilibrium turbulent flows). When (5) is thought of as the simplified 
equilibrium form of a more general pressure-strain model it becomes ambiguous as to which 
part of (5) represents the slow pressure-strain and which part of (5) represents the rapid 
pressure-strain. It can only be said definitively that eAij is the slow pressure-strain in the 
limit of relaxational turbulent flows and that KAiijkidv^/dxi is the rapid pressure-strain in 
the rapid distortion limit, 11 This ambiguity causes no problem since only the total pressure- 
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strain correlation is needed in (2) for the calculation of the Reynolds stresses. When the 
model predictions for the norm of the total pressure-strain correlation are compared with 
the DNS of Rogers et al. 8 for homogeneous shear flow a rather different picture emerges as 
can be seen in Figure 1(c). The SSG model performs as well, if not better, than the SL and 
FLT pressure-strain models. Furthermore, it was shown by Speziale, Sarkar and Gatski 11 
that none of the other models are capable of predicting the individual slow or rapid parts of 
the pressure-strain correlation for a wide range of homogeneous turbulent flows. The results 
presented in Shih and Lumley 2 are misleading in this regard. 

Comparisons such as those shown in Figure 1 are not very helpful for determining the 
predictive capabilities of a pressure-strain model in turbulent shear flows. A better test for 
gauging the performance of a model is to determine its ability to predict accurate equilib- 
rium values for the Reynolds stress anisotropies. These are the crucial physical quantities 
in homogeneous shear flow that are independent of the initial conditions and, therefore, 
repeatable. Over thirty independent test runs of homogeneous shear flow in reported phys- 
ical and numerical experiments 7-9 have yielded equilibrium values for the Reynolds stress 
anisotropies that lie within 10% of one another. On the other hand, the time evolutions of 
the individual Reynolds stresses vary by factors of 3 or 4 depending on the initial conditions 
of the test case. 

In Table 1, the equilibrium Reynolds stress anisotropies predicted by the SL, FLT and 
SSG models are compared with the experimental results for homogeneous shear flow for 
P/e = 1.5. This value of P/e is chosen since it is the average value obtained from the 
most recent experiments. The SSG model is in closest agreement with the experimental 
data whereas the SL model exhibits the largest discrepancies. The same trend is exhibited 
for other equilibrium values of P/e in the experimental range of 1.4 - 1.8 (see Abid and 
Speziale 10 for more details). 

As discussed earlier, the equilibrium Reynolds stress anisotropies for the log-layer of a 
two-dimensional turbulent boundary layer can be obtained from (9), after setting P/e = 1. 
The results obtained from the SL, FLT and SSG models are compared with experimental 
data 12,13 in Table 2. The results obtained are similar to those obtained for homogeneous 
shear flow: the SSG model is in the closest agreement with the experimental data while the 
SL model exhibits the largest deviations. In regard to the latter point, it is surprising how 
poorly the SL model performs in its prediction of the normal Reynolds stress anisotropies. 

In Figure 2, the model predictions for the time evolution of the turbulent kinetic energy 
in homogeneous shear flow are shown for two test cases: the large-eddy simulation (LES) 
of Bardina et al. 14 and the DNS of Rogers et al. 8 (run C128X). Here, the dimensionless 
turbulent kinetic energy and dimensionless time are given by K* = K [ K 0 and t” = St, 
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respectively. These results exhibit the same trend as before: the SSG model performs the 
best and the SL model has the largest deviations. The low growth rate predicted by the SL 
model arises from its underprediction of the equilibrium value of b\ 2 (see Table 1). 

This brings us to the basic question as to why the SL model performs poorly in these 
simple shear flows. As discussed in Shih et al. 1 , the SL model was, to a large extent, calibrated 
based on realizability constraints. It was recently shown by Speziale, Abid and Durbin 15 that 
the SL model yields unrealizable results in homogeneous shear flow due to an error in their 
analysis. In Figure 3, the time evolution of the invariant function F = 1 + 9/7 + 2777/ 
predicted by each model in homogeneous shear flow is shown for the anisotropic initial 
conditions: b u = -0.32, 622 = 633 = 0.16, b i2 = 623 = 613 = 0, and SK/e = 15. For 
realizable turbulence, we must have F > 0. It is clear that for these initial conditions, the 
FLT and SSG models yield realizable solutions whereas the SL model yields unrealizable 
results! Hence the primary theoretical constraint by which the SL model was formulated is 
in error. 

In conclusion, it must be said that the good performance of a model in these basic 
turbulent shear flows is no guarantee that it will perform well in more complex turbulent 
flows. However, these simple test flows do bear directly on how well a model will perform 
in equilibrium turbulent boundary layers which form a cornerstone for many engineering 
applications. A model that cannot predict these simple test cases accurately should be 
abandoned for use in more complex engineering flows. 
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APPENDIX 


The detailed form of the turbulence models considered in the paper are as follows: 
Shih & Lumley Model 

= -Cxebij + ^KS tJ + 12 a 5 K (l b ik S jk + b jk S ik 
0 



—-bkiSki^ij) + -(2 — 7a 5 )K(bikVjk + bjkSik) 
O Jo 




( Al) 


+ w zK{bubi m Sj m ~t~ bjibim S tm 2bikEkihj 
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_ 4 _ 

— 36^/ “1“ ) 
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Ci = 2 + 

y exp(-7.77/ v /.Re t ){72/ v /^ + 80.1 J»[l + 62.4(-// + 2.3///)]} 

(A 2) 


F =1 + 911 + 27 III 

(A3) 


II ~ --bijbijy III = 
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4 A 2 
Ret = q — 
9 i/e 

(A5) 


“»=s( i+ H 

(A6) 

C'i 

= 1.20, C e 2 = - + 0.49exp(— 2.83/^/fiet)[l — 0.33/n(l — 55//)] 

(A7) 

Fu, Launder & Tselepidakis Model 



n,j = — C\€bij + C*2£ 



-f — KSij + 1.2JT | 'bikSjk + bjkSik — —bkiSkifiij'j 

26 4 — 

+—K(b ik Wjk + bjkSik) + -K(b ik bkiSji 
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4 _ 14 
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C ! = -1201 1 Vf - 2s[F + 2, C 2 = 144/7\/F 
C e i = 1.44, C e2 = 1-90 

Speziale, Sarkar & Gatski Model 

n,i = -(C ie + <7^ + <7 2 e (6, fc 6 fcj - !&«&«*,•) + (<7 3 - C^)K^ {j 
+C 4 K (bikSjk + bj k Sik — -bkiSkifiij^J + C*,K{bikujjk + bjk&ik) 

Ci = 3.4, C{ = 1.80, C 2 = 4.2, C 3 = <7* = 1.30, C 4 = 1.25, C 5 = 0.40 

5 

Jib — b,j b t j 

C e i = 1.44, C c2 = 1.83. 

In these models, Sij and u7 u are, respectively, the symmetric and antisymmetric parts 
mean velocity gradient tensor dvijdxj which are given by 



- _ 1 ( dvA 
“ li ~ 2 [d^ ~d^,)' 


(-« 0 ) 


mu) 

(> 112 ) 

( 413 ) 

(414) 
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( 415 ) 

( 416 ) 
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Equilibrium 

Values 

SL Model 

FLT Model 

SSG Model 

Experimental 

Data 

bn 

0.105 

0.177 

0.214 

0.20 to 0.21 

b 12 

-0.121 

-0.157 

-0.163 

-0.14 to -0.16 

b 22 

-0.107 

-0.122 

-0.140 

-0.14 to -0.15 

b 33 

0.002 

-0.055 

-0.074 

-0.05 to -0.07 


Table 1. Equilibrium Reynolds stress anisotropies in homogeneous shear flow: Comparison i 

of the model predictions with physical and numerical experiments. 7-9 
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Equilibrium 

Values 

SL Model 

FLT Model 

SSG Model 

Experimental 

Data 

b n 

0.079 

0.141 

0.201 

0.20 to 0.23 

b 12 

-0.116 

-0.162 

-0.160 

-0.14 to -0.16 

b 22 

-0.082 

-0.099 

-0.126 

-0.13 to -0.15 

b 33 

0.003 

-0.042 

-0.075 

-0.05 to -0.10 


Table 2. Reynolds stress anisotropies in the log-layer of an equilibrium turbulent boundary 
layer: Comparison of the model predictions with experiments. 12-13 
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Figure 1. Comparison of the model predictions^ for the norm of the pressure-strain correlation 
with the direct numerical simulation of Rogers et a/. 8 for homogeneous shear flow (run 
C128U): (a) slow part, (b) rapid part and (c) total. 
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